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Abstract. Nonlinear growth of one-dimensional density structures with a frozen-in magnetic field is investigated in 
Newtonian cosmology. A mechanism of magnetic field amplification is discussed. We discuss the relation between 
the initial conditions for the velocity field and the basic time-scales determining the growth of the magnetized 
structure. 
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Q During the last decade the evidence large-scale cosmological magnetic fields has systematically grown. Fields of several 
microgauss have been measured beyond the galaxy clusters (Kim et al. 1991); the Rotation Measure (Kronberg 1994) 
confirms the existence of coherent magnetic field on Mpc scales or larger. The recent discovery of large-scale diffuse 
radio emission testifies to the presence of magnetic fields of ^ 0.1 fiG, along the 6 Mpc filament (Bagchi et al. 2002), 
with evidence for their coherent nature. 

Understanding of the magnetic field behaviour at different phases of the matter-dominated era is crucial to explain 
5J] the microgauss fields observed in high redshift objects z > 2; starting from the damped Ly a systems, through the 
5^ , distant radiogalaxies and the galaxy clusters up to the scales typical of galaxy superclusters. While the magnetic field 
of galaxies may result from dynamo effects, the fields at larger scales cannot be explained by the same mechanism. 
Here there is either no rotation, necessary for dynamo action, or the structures are dynamically too young to leave 
dynamo action enough time to operate. Magnetic fields at Mpc scales are likely to be primordial. Some pre-dynamo 
mechanisms of the primeval magnetic field amplification must be at work at least in the linear and weakly nonlinear 
regime. 

In this paper we discuss the mutual relationships between density growth and the magnetic field evolution in early 
nonlinear stages, when the wall or filamentary structure is formed. We investigate whether the growing planar density 
structure may drag and amplify the magnetic field. We employ the exact solutions for 2-D (pancake ) inhomogeneity 
evolution in the Newtonian description^ and emphasize the role of initial conditions, in particular, the large-scale 
primordial flows. The magnitude of primordial velocity fields at recombination determines the time and the growth 
rate of density fluctuation. As a consequence, it defines the duration time of the pre-dynamo and dynamo amplification 
phase. To avoid problems with a mathematical definition of the weakly nonlinear regime we work with fully nonlinear 
equations and their solutions. Although finally we refer to the regime where the density contrast A is between 1 and 
100 (which is relevant for the cosmological structures we discuss), dynamical equations are true for A > 100. 

Magnetohydrodynamic equations in the covariant notation are given in Section |2] Simplifying physically relevant 
assumptions and the resulting nonlinear perturbation equations are discussed in Section |21 Nonlinear solutions for the 
density contrast and the magnetic field enhancement are given in Section 0] Section [S] contains numerical estimations 
and graphical presentation of the magnetized pancake formation. 



^ The problem is opposite to that formulated by Wasserman (1978) and Kim et al.(1996), where the magnetic fields are 
expected to actively support the structure formation processes. 
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2. Magnetized self-gravitating fluid 

Magnetohydrodynamic equations for self-gravitating fluid with infinite conductivity (Chandrasekhar & Fermi 1953, 
Wasserman 1978, Papadopulos & Esposito 1982), expressed in a manifestly covariant way, form a dynamical system 
for the density p, the expansion rate i), and the magnetic field 

P = (1) 



i9 = r - At^kp - 2a^ + 2a;^ + -(iJ'i?,- - -H'S])'^ - -p' (2) 

3 \p 2 ' p J 

W = - (3) 

Hydrodynamic scalars ti and w, which measure the shear and the rotation, respectively, 

w = i(wifca;*'')^/^, t^ifc = ^(wi;fe - Wfc;*)- (5) 

11* stands for the fluid velocity, its divergence builds the expansion rate ■& = v\. Semicolon means the covariant derivative 
with respect to the space coordinates, dot stands for the time convective derivative, e.g. p = dp/ dt + p.^i, while (5*-' 
is the Kronecker delta; the Einstein summation convention is employed. Equations (|1I3|I are coordinate- independent. 
(In the Cartesian coordinates, covariant derivatives reduce to the partial derivatives, and tensor indices are raised 
and lowered by Kronecker delta). Equation is the continuity equation. Equation (j^J is a Newtonian analogue of 
the Raychaudhuri-Ehlers equation (Ellis 1971) — a scalar form of the equation of motion for continuously gravitating 
media with the Poisson equation included. System pi3|l has the same form in general relativity and Newtonian theory 
(Ellis 1971, compare also Tsagas & Barrow 1997), and therefore is of particular interest in the cosmological context. 

The approximation of inflnite conductivity results in a vanishing electric field. The magnetic contribution to the 
equation of motion is quadratic in iJ* (the last term in eq. Q), and therefore, can be neglected in the weak magnetic 
field limit. In this limit the system (|1I3|I splits into the autonomous system (|1I2() and the induction equation 10). Then, 
the fluid dynamics is entirely determined by gravitational forces, while the magnetic field is dragged along by fluid 
flow. 



3. Planar symmetry 

Below we consider a model based on the following assumptions: 

1. The unperturbed "background space" is an isotropic and homogeneous (Newtonian) universe. 

2. The initial perturbation state is given at random, i.e. perturbation in the density and velocity flelds l^[ini) 



Sp_ 



P(0) 



and — [^]{ini) independent quantities. The perturbation is initially small, i.e. l^[ini) <C 1 and 

(mi) 



0(jm) < 1 

3. After the recombination epoch, the infinite conductivity approximation is adequate. The electric field vanishes 
(£'* ~ 0), while the magnetic field is frozen (Chandrasekhar & Fermi 1953). The matter pressure is negligible 
(P = 0) ^ 

4. During the considered period the primordial magnetic fields and their gradients are small compared with the 
density and the density gradients, respectively {HiW <^ p, Hi-jW <C p-j). In the noncovariant approach to MHD this 
implies a small value of the ratio of Lorentz force Fl to the gravitational force Fg-"^ 

5. The perturbation is rotationless and has a planar symmetry, i.e. the velocity potential V{t, a;*) can be expressed as 
V{t, x^) = Vo(i, x^) + 5V{t, x^), where Vb(t, x*) = ^'d(Q)Xkx'' stands for unperturbed Hubble flow, and the perturbation 
5V{t, x^) is independent of two of Cartesian coordinates. Consequently, all the hydrodynamic scalars depend solely on 
time and on the only one space variable x'^ = z , the one parallel to the fluid contraction (orthogonal to the pancake 
plane). 

Under these assumptions the system (^-(El can be divided into the background dynamics 

P(o) = /3(o)^9(o) (6) 

1^(0) = -^^^(0)^ - 47rKp(o) , (7) 

^ We have Fl/Fg oc npj^p^Gpi ■f^("^i^)^' where 1 is the typical length scale and taoii — the characteristic time scale of 
collapse. Expressing then the collapse time by the quantity used in the above notation, tcoii oc jg, one obtains for the ratio: 
^ijjg)^ ^ ^ ^'-'^ relevant Mpc scales of cosmological structures. 
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where the "unperturbed" magnetic field is obviously absent, and the propagation equations for inhomegeneities 



-AnnSp - 
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(8) 
(9) 

(10) 



The perturbations Sp = p — p(o) and Si) = d — i9(o) are defined as differences between local and background values. 
They are not assumed to be small, hence nonlinear terms in 1)819(1 remain. The shear tensor for planar perturbation 



(11) 



has been employed to derive equations ij^and lTUI) . In particular the shear scalar a has been expressed by the variation 
in the expansion rate SS, i.e. a — -^S-d. The autonomous sub-system (|8l9f) can be evaluated to a single second order 

differential equation for the density contrast A — 



P(0) 



(A)- 



(1 + A) 
and solved independently. 



--i9(o)A + 47rKp(o)A(l + A) 



(12) 



4. The evolution of planar structures 

Equation (|12|1 reduces to the Jeans-Bonnor equation (Bonnor 1957, Weinberg 1972) 
2 



-■i9(o)M + 47rKp(o)M 



(13) 



under the change of the perturbation variable (compare Buchert 1989 in the context of the Zeldovich approximation) 
fi = A{l + A)-\ (14) 

This enables one to express the finite amplitude (nonlinear) perturbations as functions of the solutions to the linear 
equation (|13|1 . Equation H13|) rewritten in the conformal time rj = J a{t)~^dt takes the form 



a 2 



i\ 2 
a 



= 



(15) 



where the dimensionless scale factor 0(77), formally defined by the differential equation "doivY = 3^^^, satisfies the 
Friedman equation for a dust-filled universe 



K a' a" 

^ - ^ + 2^ =0 

and can be found in the exact form 

M.J VkA 

a U? = sm 



(16) 



(17) 



Prime in H15|l . ((16|l and all the equations below is the differentiation with respect to 77, both K and M = a^{rj)p are 
constants of motion^ . The two independent solutions of the equation (|15|l read 



M 

piiv) = -10— 



1 + - -2 + \/X77C0t 2 



Kv\ . -2 

sm 



Krj 



A'2(?7) 



-3/2 



1 + tan 



cof 



(18) 



(19) 



In general relativity K means the curvature index and is traditionally set to K = — 1,0, +1. In Newtonian cosmology K 
distinguishes between different dynamical behaviours. 
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and in the K ^ limit yield Hi{ri) = Alrj'^ and /i2('7) = M '^^^rj respectively. The arbitrary solution to is 
a linear combination 

h{t]) = ci/ii(ry) + C2H2{v) , (20) 
while the original perturbation A can be found by use of the reciprocal transform 

Although the solutions ^{rj) form the linear space (|20|l . the solutions A(77) to equation H12|l obviously do not. On the 
other hand the time evolution of the magnetic field is governed by the linear equation Hl()|l. which after investing ^ 
can be integrated exactly. The solution written covariantly reads 




^ - ^'(^)^sSmU (22) 

and can also be expressed as a function of the density contrast A(z) of the pancake structure observed at the redshift 
z (we write it in the matrix form) 

ffx(z) 

Hyi^) I = ( " ' ; . ) I 1 + A(z) I I iJ, I . (23) 

i?z(z) 

Derivation of the formula H22I23|1 does not involve the Raychaudhuri equation. Therefore, (|22I23(I is well satisfied for 
arbitrary large magnetic fields, and is independent of the p = approximation. It expresses the amplification of the 
magnetic field during the one-dimensional fiuid compression in a homogeneously expanding medium. It shows that 
the orthogonal component is systematically diluted during the universe expansion, while components parallel to 
the pancake plane increase with the local density. As a consequence, the relic magnetic field in a large-scale planar 
structure is amplified and flattened to the plane. 

All the dynamical effects of gravity, pressure or magnetic field cumulate in the fi{ri) or A(z) evolution. For the 
dust filled universe in the weak magnetic field approximation the formulae H23() . H20|) . (|18|l . H19|l and (|17() constitute 
a closed form solution for the magnetic field enhancement. 



5. Numerical estimations 



It is reasonable to view the numerical results in the K ^ universe, where the linear solution for fi{z) expressed as 
a function of the redshift z takes the form 



A 



(ini) 



l + z 



(ini) 



l + z 



l + z 



(ini) 



l + z 



+ 



e 



{ini) 



l + z 



(ini) 



l + z 



1 + Z(ini) 

l + z 



Consequently the perturbations A 



5p 



P(0) 



A(z) 



e(z) = 



H{z) 



2(A(i„j) + 36(i„i)) + 3(A(i„i) - 29(i„i)) ( -y/^^^ 



, 6 = [6] and the magnetic field H{z) read respectively 
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' {^J'-^TT^^) - 2(A(.™) + 3e(„,)) - 3(A(,„„) - 2e(,„)) (^^i^^^ 

(A(mi) + 39(i„i)) - (A(i„i) - 20(mi)) ^ y/i^^^^ 

' ^ 2(A(i„i) + 39(i„i)) - 3(A(i„i) - 28(i„i)) ^-y/i^^^^ 



z + 1 



(1 + A(z))i/, 



(ini) 



(24) 



(25) 



(26) 



(27) 



Formulae (|25I27|I show how the structure formation process depends on the initial conditions at the recombination, 
i.e. on the initial density fluctuation A(j„j) and the initial flow 0(i„i) (the velocity divergence). Setting A(-j„j) = 10~^ 
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Fig. 1. Density contrast as a function of redshift z and initial velocity flow 0(i„i). 




Fig. 2. Velocity field contrast as a function of redshift z and initial velocity flow G(i„i). 

at z = 1200, as suggested by the CMBR measurements, we draw A and 8 as functions of z and the initial flow 0(i„i) 
(see Fig. 1 and Fig. 2). 

Galactic superclusters with the density contrast A < 100 are considered as nonvirialized structures (Peacock 1999). 
Their length scales of about ^ 10 Mpc are typically expected at the moment of transition from the linear to nonlinear 
regime. The evolution of these structures relative to the initial values of the large scale inflows 6 (mi) is presented 
in Fig 3. The evolutionary paths with different 0(i„i) form horizontal lines, while the solid sloped lines mark the 
beginning (A = 1, = —0.35) and the end (A = 10, = —3) of the weakly nonlinear regime (the shaded region). The 
magnitudes of the contrasts at both characteristic moments are set to be compatible with the values obtained from the 
numerical simulation (Gramann et al. 2002). The region below the shaded region represents the linear evolution, while 
the region above — the strongly nonlinear collapse. The existence of the low z structures with A < 100, confirmed 
by the observations, when compared with the theoretical low redshift behaviour (Fig 3) favours the initial inflow 
Qi^ini) — 10~^ at the recombination epoch^. 

The initial velocity fleld determines the basic time scales for the weakly and strongly nonlinear structure formation 
respectively, and consequently, deflnes the time scales for compression and magniflcation of the seed magnetic field. 
In particular, the time of entering the strongly nonlinear regime (A « 150) is crucial for the efficiency of the dynamo 
action (e.g. Widrow 2002). The amplification degree in this process exponentially depends on the time available for 
the dynamo mechanism to operate. The pre-dynamo magnetic field compression (Lesch et al. 1995) is followed by 
the dynamo amplification, which in the case of protogalaxies results in some 10^ — lO^'' enhancement in amplitude, 
depending on the time when the field reaches the observed value. The question if the dynamo mechanism is applicable 
to galactic cluster structures is still open. Even more intriguing is the origin of the almost equipartized supercluster 
magnetic fields (~ lO"*") (Sigl et al. 1998). Fig. 3 illustrates that structures with initial inflow Q(ini) — reaches 
the weakly nonlinear regime (A ~ 1) at redshift z = 3, and leaves the weakly nonlinear phase (A ~ 100) at z = 1.3. A 
dynamo may operate then for ~ 8 x lO^yr, which (taking into account the typical dynamo time scale ~ lO^yr) results 
in the ~ 2 10^ enhancement of magnetic field. On the other hand, for the structures with the initial inflows at the 
level of 0(i„i) w 6 X 10"'^ (or 0(i„i) w 5 x 10"'^), the beginning and the end of the previrial nonlinear stage occur at 

* Note that the fluid velocities and not their divergences contribute to the Sachs-Wolfe temperature formula, therefore 
is not directly determined from the CMBR satellite measurements 
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Fig. 3. Schematic diagram of weakly nonlinear evolution of the density and velocity contrasts. For both contrasts 
the A = const and © = const sections are presented. The solid lines correspond to the initial and final states of the 
previrial nonlinear evolution of A and O. The shaded area represents all the evolutionary states between (A = 1, 
= —0.35) and (A = 10, 6 = —3). Horizontal dashed lines show the evolutionary tracks. 

z « 1.4 (z w 0.8) and z « 0.4 (z « 0), respectively, which limits the dynamo operation phase to less than 2.3 x lO^yr 

(or none). In this latter case i.e. for systems which presently attain the previrial collapse state (A ~ 100 at z ~ 0), the 
primordial magnetic field constrained at z = 10^ to magnitudes below ~ 10~^ G (Barrow et al. 1997) are pre-dynamo 
amplified by merely 2 orders of magnitude. 



We provided the nonlinear exact solutions for the density, velocity and magnetic fields for the pancake -type structures 
in the Newtonian expanding universe. The approximations of the potential velocity field and vanishing matter pressure 
have been employed. The time when the compressing flat structure enters the regime of nonlinear growth is controlled 
by the initial value of velocity field at the recombination. The structures accompanied by large hydrodynamic flows 
collapse earlier, i.e. the moment when dynamo mechanism may switch on occurs at higher redshifts, which eventually 
results in stronger magnetic field enhancement. For presently observed velocity fields, 10~^, (see, e.g. Dekel 1997) 
in supergalactic structures of 100 Mpc and A ^ 1 the initial inflows 8(i„i) ~ 10"'^ and initial magnetic fields 
~ 10~^ — 10~* Gauss are expected. The result is compatible with the simulation estimations (Gramann et 



Firm evidence of primordial magnetic flclds in structures at prc-virial stages is of particular importance, as these 
fields "remember" the initial conditions and thus set constraints on the seed magnetic fields, density and velocity fields 
at the recombination. The observational techniques become more important (Faraday rotation measurements and the 
indications coming from the propagation of cosmic radiation UHE in the Local Supcrcluster) , which potentially might 
distinguish between the large scale magnetic seed component from other magnetic fields of astrophysical origin (i.e. 
resulting from galactic dynamo, outflows from radiogalaxies etc.). The rotation measure which have the same sign 
along the Supercluster plane would suggest a coherent, relic fleld at this scale. 
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